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We study slow variation (both spatial as well as temporal) of a parameter of a system in the vicinity 
of discontinuous quantum phase transitions, in particular, a discontinuity critical point (DCP) (or 
a hrst-order critical point). We obtain the universal scaling relations of the density of defects and 
the residual energy after a temporal quench, while we also unravel the scaling of the characteristic 
length scale associated with a spatial quench of a symmetry breaking field. Considering a spin-1/2 
XXZ chain we establish how these scaling relations get modified when the DCP is located at the 
boundary of a gapless critical phase; these predictions are also confirmed numerically. 

PACS numbers: 64.60.Ht, 05.70.Jk, 75.10.Jm 


When a macroscopic system is driven from a disordered 
phase to an ordered phase by lowering the temperature 
across a phase transition, there is spontaneous generation 
of topological defects, namely, vortices, kinks and domain 
walls; this has been explained invoking upon the so called 
Kibble-Zurek (KZ) mechanism. Originally proposed in 
the context of the cosmic evolution [1, 2] and following 
the subsequent generalization to the condensed matter 
systems [3], in recent years the KZ mechanism has been 
extended to the quenches in the vicinity of a quantum 
phase transition [4-7]. In fact, the remarkable progress 
in experiments with cold atoms in optical lattices, which 
can be treated as nearly isolated systems, has paved the 
way for this recent upsurge in theoretical studies of non¬ 
equilibrium dynamics of closed quantum systems. 

A number of studies have established that, when a d- 
dimensional quantum system is driven across an isolated 
quantum critical point (QCP), by changing a parameter 
of the Hamiltonian in a linear fashion as t/tq, the density 

of defects (p) satisfies the KZ scaling Simi¬ 

larly when quenched to the QCP, the residual energy den¬ 
sity, i. e., the excess energy density over the ground state 
of the final Hamiltonian (eres) scales as 
here, v is the critical exponent of the correlation length, 
and 0 is the dynamical critical exponent associated with 
the corresponding QCP. On the contrary, when quenched 
to the gapped phase, the residual energy follows a scal¬ 
ing relation identical to that of the defect. These scaling 
formulae however get modified, sometimes drastically, in 
several interesting situations [8-13]. (see for review [Id¬ 
le]). 

The concept of the KZ scaling has also been extended 
to spatial quenches [17, 18], where a parameter is inho¬ 
mogeneous along a spatial coordinate x characterised by 
a slope Xq^ with the quantum (or classical) phase transi¬ 
tion occurring at a; = 0. The order parameter is found to 
bend smoothly over the region jxj ^ ~ ob¬ 

viously, sets a length scale characterized by quantum 
critical exponents. These predictions has been confirmed 


analytically as well as numerically in several models [17- 
19]. 

The present Letter studies temporal as well as spa¬ 
tial quenches near a quantum discontinuous (first-order) 
phase transition, which, to the best of our knowledge, has 
not received much attention in the context of quenching 
and the consequent KZ mechanism scenario; however, the 
time evolution across a discontinuous phase transition 
has been widely studied from the viewpoint of quantum 
annealing (or adiabatic quantum computation) mainly 
concentrating on the scaling of the energy gap [20, 21]. 
Although quite recently, Bonati et al. ]22] proposed a 
scaling of the length scale ( near a classical discontinuous 
transition in the presence of a temperature gradient, the 
same for a sloped symmetry breaking field is not known. 

The goal of the present Letter is to show that p and 
£res after a temporal quench across (or at) a disconti¬ 
nuity critical point (DCP) (or a first-order critical point 
[23, 24]) do indeed follow universal scaling laws, given in 
terms of the associated critical exponents derived here 
for a generic DCP. On the contrary, when a symmetry 
breaking field is spatially quenched near the DCP, a char¬ 
acteristic length scale is also shown to emerge. Even 
though similar scaling relations can be derived for clas¬ 
sical quenches, our focus will be restricted to quantum 
quenches only. 

We initiate our discussion with a phenomenological de¬ 
scription of discontinuous transitions, illustrating with 
the example of simple Ising-like systems with Z 2 sym¬ 
metry exposed to a tunable quantum fluctuation with 
its strength denoted by P. Due to Ising symmetry, dou¬ 
bly degenerate ordered ground states coexist at a low 
P in the absence of a symmetry breaking field. With 
increasing P, the order parameter characterizing the or¬ 
dered states vanishes discontinuously at a discontinuous 
transition Pc. Application of a symmetry breaking field 
h, on the other hand, lifts the degeneracy of two states 
with nonzero order parameter; as h is tuned across the 
zero-field coexisting line (see Fig. 1), the order parame¬ 
ter shows a jump from one to the other. One may now 
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FIG. 1: Schematics of two phase diagrams involving a discon¬ 
tinuous transition. 


wonder what is the fate of this coexisting line beyond Fc? 
There are indeed two possibilities: in the former, the co¬ 
existing line splits into two lines (Fig. 1(a)) when the 
order parameter shows a discontinuous jump to a finite 
value even at F > Fc with increasing the magnitude of h 
due to an interchange of a stable and a metastable states. 
The coexisting line ends at a critical point denoted by C 
in Fig. 1(a), beyond which metastable states disappear. 
This scenario can be explained within the framework of 
the phenomenological Ginzburg-Landau-Wilson theory. 
In the latter case, on the contrary, the coexisting line 
terminates at Fc (Fig. 1(b)) with no metastable state ap¬ 
pearing for F > Fc. The transition point Fc is referred 
to as a DCP and is the focus point of the present Letter. 

To enunciate the critical properties associated with 
a DCP, one proceeds with the following assumptions: 
at the DCP, the ground-state energy density e(F, h) is 
a continuous function of F and h, but its derivatives 
9e(F, 0)/9F and the order parameter m(F,0) (at h = 0) 
exhibit jumps as a function of F. We put them in a pre¬ 
cise mathematical form: 

e(Fc + 0,0) =£(Fc-0,0) (1) 

9£(Fc + 0,0) , 9£(Fc-0,0) 

aF ^ 5F ’ ^ ^ 

|m(Fc-0,0)| >m(Fc-t-0,0) =0. (3) 

On the other hand, right at F = Fc, m shows a discon¬ 
tinuous change with h at h = 0: 

m(Fc,+0) ^m(Fc,-0). (4) 

Equipped with the above assumptions, let us now pro¬ 
pose the scaling ansatz, first fixing h = 0: As usual, the 
correlation length (^) and the ground-state energy den¬ 
sity near a quantum critical point follow scaling forms. 


C |F - Fc|-^ and £(F, 0) - £(Fc, 0) « 

[25, 26], where v and z are the correlation length and 
the dynamical exponents respectively and ± stands for 
the sign of F — Fc. The first order derivative of £(F, 0) 
is written as 9£(F,0)/9F « This along 

with Eqs. (1) and (2) leads to = l/{d + z). In the 
other situation, when F = Fc and h ^ 0, the corre¬ 
sponding correlation length (^h) and the ground-state 
energy density satisfy ~ and £(Fc,h) — 

£(Fc, 0) « Recalling that m is determined 

by m(F, h) = —de(T, h)/dh, these relations together with 
Eq. (4) yields = l/{d + z). 

Henceforth, we move on to the non-equilibrium situ¬ 
ation arising due to quantum quenches. We consider a 
temporal quench of the parameter field (with h = 0) of 
the form F(t) = Fc — Fq^^; here t is the time, is a 
quench rate, and Fg is a positive constant. Consequently, 
the system is driven from the the disordered phase (for 
t < 0) to the ordered phase (for t > 0). Assuming a 
power-law scaling of the healing time r(T, 0) ^ ^(T, 0)^, 
for a time-varying parameter, we introduce the notion of 
the instantaneous healing time T{t) = ^(F(t),0)^; using 
^(F,0) ~ |F — Fcj”'^ and the quenching form of F, this 
can be written in the form r(t) ~ If the sys¬ 

tem is initially prepared in the ground state, it remains 
in the instantaneous ground state as long as the healing 
time is short i.e., the system is far away from the QCP 
at t = 0. This adiabatic time evolution breaks down, 
when t > 0 as a consequence of the diverging healing 
time. In order to derive the scaling of the healing length, 
we invoke upon the adiabatic-impulse approximation [7] 
which assumes that the the system is in the instanta¬ 
neous ground state when t < —T{t), while it freezes when 
t > —T(t). The freezing time t is then determined by 
the condition jtj = T{t). One then immediately finds 
^ ppj. f y system has an 

inhomogeneity which is characterized by a length scale 
which can be derived form the scaling relation 

6«e(F(t),0)^4/(‘'+^^^ (5) 

Noting that the defects are distributed in a d dimensional 
phase space, the scaling of the density of defects is then 
given by p ^ which yields 

If the elementary excitation in the system after a tempo¬ 
ral quench has a gapless dispersion as jkj^ with the wave 
vector k, the residual-energy density follows the scaling 
relation 


^res 


, — (^d-{-z) / (^d-\-2z) 

> 


( 7 ) 


otherwise, it satisfies the scaling given in Eq. (6). 

In a similar spirit, let us now consider a spatial quench 
of a parameter of the Hamiltonian, h{x) = ho(x/xQ), 
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keeping F fixed at Fc, where x refers to one of the d spa¬ 
tial coordinates, Xq^ is the quench rate or the scale of in¬ 
homogeneity, and ho is a positive constant. Note that the 
system is critical when the field vanishes, i.e., at around 
a; = 0. Defining the local correlation length[\?>\ ^(x) = 
^(Fc,/i(x)) and using the relation ^{Tc,h) ~ \h\~‘'^ 
(^'H = z)), one readily finds ^(x) ~ (|x|/xq)“'^“. 

In spite of the inhomogeneity, the system should main¬ 
tain the local equilibrium for x sufficiently far from x = 0 
since the local correlation length becomes is much smaller 
than XQ. On approaching the critical situation, however, 
the local correlation grows and hence the approximation 
of the local equilibrium breaks down leading to a finite 
correlation length near x = 0. To make an estimate 
of this finite correlation length, one again invokes upon 
the pulse-impulse approximation: the system maintains 
its local equilibrium for x satisfying |x| ^ while 

for |x| <C ^(x) the system freezes with the correlation 
length The boundary x of these two regimes is ap¬ 
proximately determined by the condition f{x) = |x| and 
f is given by = £,{x) which results in the scaling rela¬ 
tion |x| ~ = XgOne therefore arrives 

at the scaling relation of the length scale 

(8) 

Equations (5), (7) and (8) constitute the main result 
of our paper regarding the quenching in the vicinity of a 
generic DCP as depicted in Fig 1(b). We generalize these 
relations below by exploring the dynamics of a XXZ chain 
in a magnetic field described by the Hamiltonian 

L-l L 

+ ^«+i) - hiaf, (9) 

i=l i=l 

where erf’s (o? = x,y,z) are the Pauli operators and L 
is the system size. We first set hi = 0,Vz, when the 
model exhibits a discontinuous quantum phase transi¬ 
tion at A = — 1. The ground state for A < —1 possesses 
the complete ferromagnetic order, while it is a disordered 
state for —1 < A < 1. Starting from the the ground 
state of the initial Hamiltonian with A = 0, we con¬ 
sider a slow ramp to A = — 1 following a linear protocol 
A(t) = —t/tq. Since the Hamiltonian commutes with 
'^tot = the initial state lies in the erfot = 0 

sector, the dynamics is restricted within this manifold 
only. Even then the transition occurring at A = Ac = — 1 
indeed can be viewed as a DCP as we have established 
rigorously in the supplemental section. However, this 
DCP occurs at the end of a gapless critical line, which 
results in modification of associated critical exponents 
and consequently, of the scaling relations of p and £res- 

The low-energy states of the XXZ chain with — 1 < 
A < 1 are described by the Tomonaga-Luttinger liquid 
with gapless excitations [27, 28]; this gaplessness, rules 
out the possibility of defining the exponent v using the 


conventional scaling of the energy gap AE ~ (A — Ac)^*^ 
[25]. However, the dispersion relation at Ac is quadratic, 
i.e., energy gap AE ~ L~^ implying z = 2. 

To estimate the exponent v, let us recall that the 
ground-state energy density in the = 0 manifold be¬ 
haves as £(A) — £(Ac) = —(Ac — A) for A < —1 and 
« (A — Ac)^/^ for A > —1 [29]. Comparing with the 
conventional scaling of the same, namely, (A — 
along with d = l,z = 2, one finds that v = \/2 which 
clearly differs from the predicted value v = l/{d + z). 
This mismatch stems from the existence of the gapless 
line for A > Ac characterized by gapless excitations with 
dispersion lfcl^+, with 2 + = 1 in the present case. 

To settle this anomaly, we propose the following modi¬ 
fied scaling forms of the spectrum AE ~ jA — -I- 
and the ground-state energy density £(A) — £(Ac) ~ 
jA — Acl'^^(l -I- jA — Acl“'^)“‘’*, where z± characterize 
the excitations on either sides of the DCP respectively; 
z+ = 1 and Z- = 0 (on the gapped side), in the present 
case. Similarly, A± are the critical amplitudes (of f) on 
the gapless and the gapped sides. The case = 0 im¬ 
plies a perfectly ordered ground state without any fluc¬ 
tuation; this is indeed true for the XXZ chain which can 
be established using the same scaling of the ground-state 
energy density and the gap. 

In order to derive the critical exponents of the DCP 
associated with a gapless phase, it is convenient to in¬ 
troduce a quantity (/>(A) = (£(A) — £(Ac))lA — Acl“‘^^± 
with for A < Ac and z+ for A > Ac, and enforce 
the modified conditions for the DCP in terms of 4>{A): 
(i) follows the same scaling law for A < Ac and 

A > Ac, (ii) (/>(A) is continuous at A = Ac, and (hi) 
is discontinuous at A = Ac (see Supplemen¬ 
tal Material for details). With these conditions together 
with the assumption that the gap closes on approaching 
the DCP as jA — Ad*^^, one immediately finds: when 
0 = A- < H+, Z- = d — z+, 6± = ^{z — z±), and 
n = l/{d-\-z—Zj^) = l/{z—z-), while if A± > 0, z- = z+, 
0± = n(z — z±), and ly = l/(d + z — z+). Remarkably, 
this reduces to ^ = 1/2 for the DCP in XXZ model where 
A- = 0, A+ > 0 and z+ = 1, and n = l/{d + z) as pre¬ 
dicted earlier when A± > 0 and = 0. 

The modified scaling relations P and £res, can then 
be readily arrived at. We find 


7 ,l/(d+2z-z+) 

' .-d/{d+2z-zj^) 

P 


( 10 ) 

( 11 ) 


and 


ere. ~ (12) 

respectively. For the DCP at Ac, 0 + = 1/2, z = 2 and 
z+ = d = 1, one gets £res ~ We recover Eqs. (5) 

and (7) when z_|_ = 0 which also implies 6± = vz. 

This scaling relation is verified by numerically study¬ 
ing the time evolution using the time-dependent density 
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FIG. 2: Residual energy per spin of the XXZ chain after a 
temporal quench A(t) = —t/tq with time t from 0 to 1 and 
rate 1 /tq in the absence of a magnetic field. Results of DMRG 
for sizes L = 100 and 200 and a fitting line are shown. Fitting 
is done by assuming a function atg*’ with the data of L = 100 
and tg’s from 10 to 100 and yields a — 0.0276 and b — 0.779. 
Inset shows the scaling of Sres for a sudden quench of small 
magnitude <5 starting from the gapless phase and ending at 
the DCP. 

matrix renormalization group (DMRG) for finite systems 
with the open boundary condition. We kept m = 50 
states for each left and right blocks in DMRG and chose 
the Trotter slice At = 0.025 [32] with the initial state 
prepared by the static DMRG. Figure 2 shows the results 
for sizes L = 100 and 200, confirming that e^s indeed ex¬ 
hibits a power-law scaling with an exponent —0.779; this 
is consistent with the value —3/4 [33]. 

If the system is quenched across the QGP to the gapped 
ferromagnetic phase, the scaling of Sres (which is identi¬ 
cal to that of p for such a quench) satisfies eres ~ 
this is in perfect agreement with the result of Pellegrini 
et al. [31] for a similar linear quench of A. For a sudden 
quench of small magnitude {6) from the the gapless phase 
to the DGP, we find £res (Inset Fig. 2). We empha¬ 

size that this is in perfect agreement with the adiabatic 
perturbation theoretic prediction [34] ej-es with 

1 / = 1/2 and z = 2; this implies that dynamical exponent 
associated with the DGP dictate the scaling of £res- 

We next move to the spatial quench with A fixed to 
— 1 and estimate izh- In the homogeneous system with 
hi = the ground state with h = 0 as well as h ^ 0 
is the fully polarized state and the energy density of the 
ground state scales as £(/i)—e(0) = —\h\. The energy gap, 
on the other hand, scales with h as AE ^ \h\ ~ 

[35]; with z = 2, this immediately implies izr = 1/2. We 
emphasize here that this value of t'n differs from l/{d+z) 
predicted before. This is due to the fact that the ground 
state is completely spin polarized without fluctuations 
resulting in a that is independent of h and always 
~ 1. Consequently, the the ground-state energy den- 


FIG. 3: Local magnetization obtained by the static DMRG of 
the XXZ chain with A = — 1 and with a sloped magnetic field 
h = (i — 1/2 — Ll2)lxq. The abscissa indicates the position 
scaled by Xq . Results for several xq collapse into a single 
curve thereby confirming the predicted scaling. 

sity does no longer scale as ^ rather scales 

as using this scaling one arrives at zzh = 1/z and 

Xq = Xq ^ with z = 2. 

Let us now apply a sloped field hi = (i—1/2 — L/2)/xq 
choosing an even L and probe the local magnetization, 
rrii = (erf) of the ground state. The scaling theory pre¬ 
sented above predicts that the characteristic length scale 
near i = {L + l)/2 is given by |s ~ In fact, one 

then expects the scaling rrii ^ J- I — \j^\ ■, where J- is 

the scaling function; remarkably this prediction perfectly 
matches with the numerical results obtained using static 
DMRG as presented in Fig. 3. 

In conclusion, we provided universal scaling relations 
for quenches in the vicinity of a quantum DGP for both 
spatial and temporal quenches. Furthermore, we have 
established that indeed there exists a DGP in the phase 
diagram of a spin-1/2 XXZ chain. The critical expo¬ 
nents and hence the scaling relations following quenches 
get non-trivially modified since this DGP occurs at the 
terminating point of a gapless critical line. Most impor¬ 
tantly, we posit generic scaling relations which perfectly 
reduce to both the situations (Ising and XXZ) consid¬ 
ered here. We emphasize that the scaling predictions for 
the XXZ chain are verified numerically. An experimental 
demonstration of our results is expected in a cold atomic 
system in an optical lattice or a trapped-ion system. 
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Here, we shall establish that the quantum critical point appearing at the boundary of the gapless critical line and 
the gapped ferromagnetic phase of the spin-1/2 XXZ chain (Hamiltonian (9) of the main text) with hi = 0,Vz , 


L-l 


^ = E 


'i+l 


+ 


i+l 


Aaf 


D+l) 


(SI) 


at A = Ac = —1 is indeed a discontinuity critical point (DCP) [SI, 2] and show how the associated critical exponents 
get altered in this continuous symmetry model where the DCP is located at the boundary of a gapless critical line; in 
the process, we shall also touch upon the scaling relations and critical exponents associated with a generic DCP. The 
complete phase diagram of the chain is shown in Fig. SI [S3]. 

To this end, we recall that the ground-state energy density £(A) in the = 0 manifold is given exactly by [S4] 


£(A) 


A 


cos /i — sm fj, 


dx 


_ao cosh7ra;(cosh2y;ix — cos^) 


(for A < —1) 

(for -1 < A < 1) 


(S2) 




where pi is defined through the relation cos/i = A. When 0 < A -|- 1 <C 1, making an expansion cos pi 
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FIG. SI: The phase diagram of Hamiltonian (SI); the gapless region lying between —1 < A < 1 is the Tomonaga-Luttinger 
liquid (TLL) phase. The transition between the gapless phase and the antiferromagnetically (AFM) ordered phase at A = 
1 is Berezinski-Kosterlitz-Thouless (BKT) transition point while the critical point separating the gapless phase from the 
ferromagnetic (FM) phase at A = —1 is a discontinuity critical point (DCP). 


(which implies « A + 1) and retaining lowest order terms in Taylor expansion in 6, we arrive at 

e(A) « - 1 - 

= -1 + ^(A + 1)3/2 + 0 ((A + 1)2) . 

Consequently, one gets 

{ A — Ac (for A < —1) 

9 /9 ’ 

^(A-Ac) 3/2 (for-l<A<l) 

where we have used Ac = —1. This immediately reveals that the derivative of £(A) —e(Ac) is discontinuous at A = Ac, 
and establishes that this point is a DCP. 

Now, in order to settle the value of jz, we compare e(A) — e{Ac) ~ (A — Ac)^/^ for A > Ac with the conventional 
scaling of the ground state energy density[S5, 6], (A — Ac)/^^''"^^'^ . This comparison leads us to the relation {d + z)v = 
3/2. Note here z = 2, since the spin wave accounts for the excitation at A = Ac even in the tr/cit = 0 manifold and 
the spectrum is quadratic in the momentum. Therefore, taking d = 1 into account, one finds v = 1/2. Surprisingly, 
this is inconsistent with the prediction 1/{d -\- z) = \/?>. As we shall elaborate below, this anomaly emerges from the 
existence of a critical line for A > Ac. 

To modify the scaling theory accordingly, we assume that the spectrum follows the form | A — Ac|®±/c*± + where 
(0_,2_) and (0+,z+) represent exponents on the lower-A and the higher-A sides of the DCP, respectively. We here 
demand 0 < z± < z. Clearly, one can readily retrieve the quadratic scaling of the gap right at the DCP by setting 
A = Ac. Comparison of the spectrum with the conventional scaling form |A — Ad^^ yields 

d± = jy{z-z±). (S4) 

Next, we assume that the correlation length has the form ^(A) ^ 1 + | A — Acl”"^, where ± corresponds to the sign 

of A — Ac. It is to be noted that ^(A) ~ 1 when = 0 (or A+ = 0); this situation arises when the fluctuation of 
the order parameter is entirely absent which indeed happens, as we shall show below, in the ground state of the XXZ 
chain for A < Ac. Having assumed the the scaling form of the spectrum and the correlation length, the ground-state 
energy density can be written as 

£(A)-£(Ac)~ |A-Ac|®±+^^±e(A)-‘^~ |A-Ac|®±+‘^^±(l-b4l±|A-A]-'^)-'^. (S5) 


We emphasize that, in general, this quantity does not scale with an identical exponent on the lower- and higher-A 
sides. To resolve this asymmetry, let us introduce the new quantity: 


ip{A) 


(£(A) - £(Ac)|A - Acl'"^-) for A < Ac 
(£(A) - £(Ac)|A - Ac|''^+) for A > Ac ’ 


(S6) 


and demand modified conditions of the DCP as follow: (i) this quantity follows the same scaling relation on both 
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sides of the critical point and that (ii) it is continuous while (iii) its derivative is discontinuous at the critical point: 


(i) 

(p{Ac — e) ~ + e) ~ e “ with an a > 0 , 

(S7) 

(ii) 

(p-{Ac - 0 ) = V3+(Ac -b 0 ), 

(S 8 ) 

(iii) 

dip-{Ac - 0 ) ^ dip+{Ac + 0 ) 

dA ^ dA ' 

(S9) 


We hereafter discuss the cases with A- = 0 and A- ^ 0 separately, with A+ 7 ^ 0 in both the situations. 
We first consider the situation where A- = 0. Equations (S5) and (S 6 ) then yield 


V3(A) 


(Ac - A)^- for A < Ac 
(A - Ac)®++‘^'^ for A > Ac ’ 


(SIO) 


Using Eq. (S4) and the condition (i) on V 3 (A), one obtains {z — Z-)^ = {d + z — which immediately reduces 

to z_ = z+ — d. The condition (ii) along with (i) result in constraint relations, Z- < z and z+ < d + z (though 
these are always satisfied under the restriction we have imposed on z±). Finally, by the condition (iii), one obtains 
V = l/{d + z — z+) = l/(z — Z-). It is noteworthy that we are now left with only two independent exponents out 
of the six exponents 2;, v, z±, and involved in the scaling. For instance, if z and z_|_ are known for a given d, the 
other exponents are automatically fixed. Now let us return to the XXZ chain. The ground-state energy density and 
the energy gap of the chain scale identically as |A ~ Ad for A < Ac. This scaling together with the the fact that the 
spectrum is gapped in this region imply A- = 0 and establish that the XXZ chain indeed corresponds to the present 
situation. Applying the value d = z+ = 1 and z = 2 of the XXZ chain to the generic relations, one obtains ly = 1/2, 
z_ = 0, = 1, and 0+ = 1/2, which are completely consistent with the property of the XXZ chain. 

Next, we consider the case with A_ 0. Equations (S5) and (S6) now reduce to 


(/?(A) 


(Ac-Af-+‘^‘' for A < Ac 
(A - Ac)®++'^^ for A > Ac ■ 


(Sll) 


Following the same set of arguments as in the previous situation results in z_ = z+ and ly = l/(d-|-z — z+) = 
l/(d -I- z — z_). If one sets z+ = 0, one obtains ly = l/(d + z). Hence the correlation length exponent of a generic 
isolated DCP, mentioned in the main text, is retrieved. 

To summarize, we have provided the scaling theory for a generic DCP based on the conditions mentioned above. 
When a DCP characterized by the dynamical exponent z (> 0) is located between two critical phases with different 
dynamical exponents z_ and z+ (assuming 0 < z_ < z+ < z without loss of generality), the exponents z_ and z+ 
are constrained by z_|_ — z_ = d and the correlation length exponent is given by ly = l/(d + z — z+) = l/(z — z_). 
Moreover, the critical amplitude of the correlation length satisfies the condition A_/A_(. = 0. Note that these rules 
are valid even when one of the two phases is gapped (z_ = 0). On the other hand, when a DCP with z > 0 is located 
between two critical phases with an identical exponent z_ = z+ = z' < z or two gapped phases (z_ = z+ = z' = 0 ), 
one finds ly = l/(d + z — z') and the critical amplitudes of the correlation length are finite on the both sides of the 
DCP. 
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